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ABSTRACT 

The ultraluminous broad absorption line quasar APM08279+5255 is one of the most 
luminous systems known. Here, we present an analysis of its nuclear CO(l-O) emis¬ 
sion. Its extended distribution suggests that the gravitational lens in this system is 
highly elliptical, probably a highly inclined disk. The quasar core, however, lies in the 
vicinity of naked cusp, indicating that APM08279+5255 is truly the only odd-image 
gravitational lens. This source is the second system for which the gravitational lens 
can be used to study structure on sub-kpc scales in the molecular gas associated with 
the AGN host galaxy. The observations and lens model require CO distributed on a 
scale of ~ 400 pc. Using this scale, we find that the molecular gas mass makes a signif¬ 
icant, and perhaps dominant, contribution to the total mass within a couple hundred 
parsecs of the nucleus of APM08279+5255. 
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1 INTRODUCTION 

Identified serendipitously in a search for high latitude car¬ 
bon stars, the 2 = 3.9 broad absorption line quasar 
APM08279+5255 is coincident with an IRAS source with 
a flux of 0.95Jy at 100/rm (Irwin et al. 1998). Observations 
with SCUBA reveal that APM08279+5255 possesses a sig¬ 
nificant submillimetre flux of 75mJy at 850/rm (Lewis et al. 
1998), implying a bolometric luminosity of ~ 5 x 10 15 Lq. 
Imaging reveals that APM08279+5255 is not point-like, but 
rather is extended over a fraction of an arcsecond with a 
structure indicative of gravitational lensing (Irwin et al. 
1998). The composite nature of APM08279+5255 was con¬ 
firmed in adaptive optics (AO) images obtained by Ledoux 
et al. (1998), with the system appearing as a pair of point¬ 
like images separated by 0.4 arcsec. Observations with NIC- 
MOS on the Hubble Space Telescope (Ibata et al. 1999) 
and AO images taken with the Keck telescope (Egami et al. 
2000) also uncovered a fainter third component between the 
brighter two. Gravitational lens models derived from these 
observations suggest that the quasar continuum source has 
been magnified by ~ 90. 

Using IRAM, Downes et al. (1999) detected emission 
in CO (4-3) and CO (9-8), revealing the presence of warm 
circumnuclear gas in APM08279+5255. Papadopoulos et al. 
(2001) were able to search for CO(l-O) and CO(2-l) in this 


system using the Very Large Array (VLA). Both were clearly 
detected associated with the quasar nucleus, as well as a 
more extended component located several arcsecs from the 
quasar images. Using locally established values of the CO- 
to-H 2 ratio, this lone cloud represents ~ 10 11 Mq of cold 
and/or subthermally excited gas. 

In this paper, we present an analysis of nuclear CO(l-O) 
emission in APM08279+5255 using VLA [] at high spatial 
resolution (0.3 arcsec). The CO appears as a partial ring of 
~0.6 arcsec diameter. These data suggest a total revision in 
the gravitational lens model for this source, with the new 
model involving a ‘naked cusp’, which naturally accounts 
for the observed odd-number of images. They also imply 
that the nuclear CO must be spatially extended on a scale 
of at least 400 pc, making this the second source in which 
gravitational lensing can be used as a ‘telescope’ to explore 
sub-kpc scale structure of molecular gas in the AGN host 
galaxy. 


1 The VLA is operated by the National Radio Astronomy Ob¬ 
servatory, which is a facility of the National Science Foundation, 
operated under cooperative agreement by Associated Universities, 
Inc. 
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Figure 1. The 23Ghz continuum (left) and CO(l-O) line emission (right) in APM08279+5255. The contours are at -1 (dashed),0,1... 
Xcr. The image has been CLEANed and the ellipse in the upper corner of the left-hand panel represents the CLEAN beam. The quasar 
A image, as determined from NICMOS imaging, lies at 08 31 41.64, 52 45 17.5, offset from the radio emission (Ibata et al. 1999). This 
is probably due to astrometric errors in the VLA and HST reference frames. 


2 OBSERVATIONS AND RESULTS 

Observations of the CO (1-0) emission from 
APM08279+5255 were made in March and April 2001 
using the VLA in the B (10 km) configuration. A total 
of 20 hours were spent on the source. Due to limitations 
with the VLA correlator, we chose to observe in continuum 
mode using two 50 MHz bandwidth IFs, each with two 
polarizations. One IF was centered on the redshifted 
CO(l-O) line, corresponding to a frequency of 23.465GHz, 
while the second IF was tuned away from the emission 
line to a frequency of 23.365GHz. This observing set-up 
maximizes sensitivity since the effective bandwidth (45 
MHz, 575 km/s) closely matches the obseved CO linewidth 
(Downes et al. 1999), but sacrifices velocity information 
(Carilli, Menten, & Yun 1999). The source 3C286 was used 
for absolute gain calibration. The rms noise on the final 
image is 35/rJy beam -1 with an effective spatial resolution 
of FWHM=0.39" x 0.28" with a major axis position angle 
of —70°. Fast switching phase calibration was employed 
with a calibration duty cycle of 130 seconds (Carilli & 
Holdaway 1999). The phase stability for all the observations 
was excellent, such that the array was easily phase coherent 
over the calibration cycle time. This was demonstrated by 
imaging the phase calibrator (0824+558) with a similar 
calibration cycle time, and by imaging the radio continuum 
emission from APM08279+5255 itself. 

The images of the line and continuum emission from 
APM08279+5255 are shown in Figure 1. The CO line im¬ 
age was generated by subtracting from the visibility data a 
Clean-Component model of the continuum emission made 
from the off-line data. The peak surface brightness for the 
continuum emission is 0.24 ± 0.035 mjy beam -1 , with a 
total flux density 0.41 ± 0.07 mjy. The corresponding num¬ 
bers for the line emission are 0.183 ± 0.035 mjy beam -1 
and 0.39 ± 0.09 mjy. The latter corresponds to a velocity- 
integrated line flux of 0.22 + 0.05 Jy km/sec, consistent with 
0.15 ± 0.045 Jy km/sec deduced for the inner ~ larcsec by 
Papadopoulos et al. (2001) 

The continuum emission is extended north-south, with 
a position angle and spatial extent as expected based on the 
three component optical gravitational lens. An attempt was 
made to decompose the continuum map into three point¬ 


like images, centred upon the positions derived from the 
analysis of HST images (Ibata et al. 1999). Given the im¬ 
age resolution, it was not possible to cleanly separate the 
north-most images (A+C) which have a separation of ^0.15 
arcsecs. The relative flux of (A+C)/B~1.5, quite similar to 
the optical ratio presented by Ibata et al. (2001). The line 
emission, on the other hand, shows distinct curvature away 
from the continuum position angle, and at the 3<r surface 
brightness level it appears as an almost complete ring with 
a diameter of about 0.6 arcsec. While the map of the line 
emission is of low signal to noise, the ring-like structure is 
apparent in each of several days worth of observations. Also, 
its size is in agreement with that of the CO (4-3) and CO (9- 
8) emission (Downes et al. 1999) and the CO(2-l) emission 
(Papadopoulos et al. 2001). Hence, these observations reveal 
the gross features of the CO (1-0) emission, although more 
observations are required to uncover the finer details. 

The observed CO (1-0) brightness temperature of the 
ring (averaged over ~575 km/sec) is 1.4+0.8 K, which cor¬ 
responds to an emitted brightness temperature of T), ~7 
K at z~3.9. A lower limit for the magnification factor of 
the CO(l-O) emission can be derived assuming the warm 
gas emitting the high-J CO lines (Downes et al. 1999) to 
be also emitting the J=l-0. This gas phase is optically thick 
with Tkin ~ 200A', in agreement with the inferred dust tem¬ 
peratures (Lewis et al. 1998), which then yields a velocity- 
averaged filling factor of / ~ 7/200 = 0.035. Since differ¬ 
ential lensing will “boost” a compact warm region at the 
expense of a more extended and possibly sub-thermally ex¬ 
cited gas phase that emits the CO J=l-0, the true value of 
f will be larger. 


3 GRAVITATIONAL LENSING 

The lensing galaxy has yet to be identified in 
APM08279+5255. Hence, only the quasar positions 
and magnitudes are available to constrain any lensing 
model. Other than APM08279+5255, all other lens systems 
possess an even-number of images. Theoretically, any 
non-singular mass distribution should produce an odd 
number of lensed images (Burke 1981), and the ubiquity 
of even numbers of images has been used to limit the 
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core radius in lensing systems, as small cores result in the 
demagnification of one of the images. Due to the brightness 
of the central source, however, models for APM08279+5255 
have required the opposite, a very circular model with a 
large core / shallow cusp (Lewis et al. 1999; Egami et al. 
2000; Munoz et al. 2001). We further explore the lens model 
of APM08279+5255 in light of the observations presented 
here. 

The concordance between the optical images of 
APM08279+5255 and the radio continuum at 3.6cm (Ibata 
et al. 1999) and 23GHz demonstrates they come from coin¬ 
cident regions with a similar scale size. Associated with the 
active core of APM08279+5255, these regions are smaller 
than the gravitational lensing caustics. The quite different 
morphology displayed in the CO image, however, indicates 
that this emission arises in a larger region, distinct from 
the continuum radiation. With its extended nature, the CO 
emitting region can lie under different parts of the caustic 
network. The CO source, however, can not be significantly 
more extended than the caustic network, as the resulting 
image would not show the ring structure seen in Figure OJ 
Therefore, for a given model, the size of the ring image pro¬ 
vides a probe of the scale size of the emitting region. 

To form ring-like images, the source must be extended 
and cover a substantial fraction of the inner caustic structure 
of the lens; the resulting image is appears as a ring, follow¬ 
ing the outer critical line in the image plane (see Kochanek, 
Keeton & McLeod 2001). Examining the critical line struc¬ 
ture in all previously published models for APM08279+5255 
(Ibata et al. 1999; Egami et al. 2000; Munoz et al 2001), the 
resulting structure for the Einstein ring of a source centred 
upon the quasar nucleus should be quite circular, with a 
radius of ~ 0.2 arcsecs, passing through the two brighter 
quasar images; such images for extended source can be seen 
in the models of Egami et al. (2000). This is quite different 
to the CO structure displayed in Figure p], with the CO emis¬ 
sion clearly showing a roughly east-west extension, with the 
hole in the ring occurring ~ 0.5 arcsecs from image A. Using 
the model derived from the HST data (Ibata et al. 1999), 
we explored the image configurations a range of sizes and 
shapes, centred upon the quasar source, for the CO emit¬ 
ting region. None reproduced the observed image structure 
and we reject this previous model. 

Naked cusps occur in highly elliptical systems, such as 
flattened disks, when the inner diamond caustic extends out¬ 
side the elliptical caustic. A source inside this extension 
produces three roughly colinear images of similar bright¬ 
ness (Mailer, Flores & Primack 1997; Bartelmann & Loeb 
1998; Keeton & Kochanek 1998; Moller & Blain 1998; Blain, 
Moller & Mailer 1999; Bartelmann 2000). While gravita¬ 
tional lensing by spiral galaxies has been observed (e.g. 
B 1600+434; Koopmans, Bruyn & Jackson 1998) no ob¬ 
served lensed quasar system has so far been associated with 
a naked cusp. 

To examine the possibility that the observed image con¬ 
figuration is due to gravitational lensing by a naked cusp, we 
constructed a simple mass model. Following Batelmann and 
Loeb (1998), this consists of a truncated flattened disk in 
a spherical halo. APM08279+5255 has brightened consider¬ 
ably since its disc overy (Lewis, Robb fe Ibata 1999; E. Ofe k 
priv. comm., see ittp: //wise-obs . tau. ac . il/~eran/Ll ), 
potentially due to the effects of gravitational microlensing. 


Comparing the images obtained in 1998 (Ibata et al. 1999; 
Egami et al. 2000) and those obtained in 1999 (Munoz et 
al. 2001), the relative image brightness have changed ap¬ 
preciably, with image B changing from being 78% of image 
A to only 50% in the latter epoch; such behaviour is ex¬ 
tremely suggestive of the action of gravitational microlens¬ 
ing, although longer term monitoring is required to confirm 
this. Hence, the relative image brightnesses cannot be used 
to constrain any mass model. Given the sparsity of con¬ 
straints, we choose to find a model that can recover the 
image positions, while providing a reasonable description of 
the observed CO emission. As can be imagined, the param¬ 
eter space is large, so a range of models that reproduced 
the quasar positions were chosen and then modeling of the 
CO emission was undertaken by-eye. With this, therefore, 
we do not claim that the model presented here is unique, 
only consistent with the general form of data. 

In our chosen model, the disk is highly inclined, present¬ 
ing a projected axis ratio of 0.25, and the disk truncated at 
8/i _1 kpcrl and possess a core radius of 0.065h^ 1 kpc. The ro¬ 
tation velocity of the disk is 200kms~ 1 . Figure ^ presents the 
source and image plane for this model, with the elliptical and 
diamond caustic, and corresponding critical lines apparent. 
In this model, the quasar images are not substantial mag¬ 
nified, with a total magnification of ~ 7, with the intrinsic 
source of APM08279+5255 being correspondingly luminous, 
Lbol ~ 7x 10 14 I/0. While extreme, this value is not necessar¬ 
ily outrageous as the unlensed quasar HD1946+7658 possess 
an intrinsic luminosity of ~ 4 x 10 14 L@ (Hagan et al. 1992), 
and APM08279+5255 may be a member of this very lumi¬ 
nous class of quasars. It must be conceded, however, that 
the non-uniqueness of the lens model translates into uncer¬ 
tainty in the model magnification and a true determination 
of the intrinsic properties of APM08279+5255 require mod¬ 
els derived from better observational constraints. 

The CO source is taken to be have an circular sur¬ 
face brightness distribution centred upon the quasar posi¬ 
tion. One important aspect of the results presented herein 
is that APM08279+5255 becomes the second system for 
which the gravitational lens can be used to study struc¬ 
ture on sub-kpc scales in the molecular gas associated with 
the AGN host galaxy, the first system being the Clover Leaf 
quasar, H1413+117 (Yun et al. 1998; Kneib, Alloin, & Pcllo, 
R. 1998). For APM08279+5255, the observations and lens 
model require the CO to be distributed on a scale covering 
a substantial fraction of the caustics in the image plane, but 
not too large to lose the ring structure. The lower limit to 
the CO source size based on the modeling is ~ 400/i -1 pc, 
while a rough upper limit is ~1 kpc. With this model, the 
CO(l-O) has been magnified by a factor of ~ 2.5 — 3. Like 
the Clover Leaf, we find that the spatial extent and mass 
of the molecular gas in APM08279+5255 are comparable to 
those seen in nearby nuclear starburst galaxies (Sanders and 
Mirabel 1996; Downes and Solomon 1998). 

Downes et al. (1999) determined a CO source size of 
~ (80 — 135)ft _1 pc for the estimated magnification factors 
of ~ 20 — 7. This size is much smaller than the one calculated 
above, while their magnification factors are larger. Their 
analysis is based upon modeling of CO emission in the grav- 


2 Q = 1 and A = 0 is assumed throughout 
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Figure 2. The source plane (left) and image plane (right) for the gravitational lens model for APM08279+5255. The diamond and 
elliptical caustics, and their corresponding critical lines, are apparent in the panels. The contour lines represent the CO source and image 
and the filled circles represent the positions and brightnesses of the quasar images as determined from the HST images (Ibata et al. 
2000). Note, however, the model presented in the paper results in brightness ratios of C ~ B ~ 0.75/1, whereas the observed ratios are 
C ~ 0.2 A and B ~ 0.8A 


itationally lensed ultraluminous galaxy IRAS F10214+4724 
(Downes, Solomon & Radford 1995), whose image is clearly 
an extended arc-like feature which possesses an essentially 
linear magnification. Such a simple model is probably a poor 
representation of the lensing in APM08279+5255. Addition¬ 
ally, Downes et al. (1999) assumed that the velocity filling 
factor is unity, substantially larger than the value derived 
in Section as the intrinsic source radius in their model 
scales inversely with this value and the magnification fac¬ 
tor is proportional to it. For f~0.35, a value well within our 
estimated range, the Downes et al. (1999) model yields an 
upper limit for the intrinsic source size and a lower limit for 
the magnification factor that our similar to ours. 

Accounting for the influence of gravitational lensing, 
the velocity-integrated CO (1-0) flux density implies that 
the nuclear content of molecular gas in APM08279+5255 
is ~ 10 10 h~ 2 Mq, assuming a CO-to-Fh conversion fac¬ 
tor of ~ 1(M@ km sec -1 pc 2 )^ 1 which is typical for star- 
burst /kinematically violent, UV-intense environments of 
gas-rich, IR-ultraluminous systems (Downes & Solomon 
1998). Assuming the CO is in a rotating disk, the dynam¬ 
ical mass can be calculated from the radius of ~ 500 pc 
set by the lens modeling, and using a rotational velocity of 
350 km sec -1 set by the observed of line velocity HWHM 
= 250 km sec^ 1 and assuming a disk inclination angle of 
45° (Downes et al. 1999). The implied dynamical mass is 
1.5 x 10 10 Mg within ~500 pc of the nucleus, consistent with 
the value derived from the CO flux. Hence it appears likely 
that the molecular gas mass makes a significant, and per¬ 
haps dominant, contribution to the total mass within a few 
hundred parsecs of the nucleus in APM08279+5255, unless 
the nuclear CO disk is close to face-on. A similar conclusion 
has been reached for most nearby nuclear starburst galaxies 
(Downes and Solomon 1998). 


4 CONCLUSIONS 

This paper has presented resolved images of nuclear CO(l- 
0) emission in the gravitationally lensed BAL quasar 
APM08279+5255. While the continuum emission is found to 


be well aligned with the optical quasar images, the CO (1-0) 
is more extended, with a broken ring-like appearance. Such a 
structure is consistent with the action of gravitational lens¬ 
ing, with the continuum emission occurring on the scale of 
the quasar core, while the CO (1-0) arises from a larger re¬ 
gion and is differentially magnified. The three-image nature 
of APM08279+5255 has posed a problem for lens modeling, 
as an extremely large, flat core is required to produce the 
central image. Such three image configurations are a nature 
consequence of gravitational lensing by a flattened potential 
which can produce naked cusps. Modeling of the CO (1-0) 
emission supports this hypothesis, although a deficit in con¬ 
straints implies that the model is not unique. An immediate 
prediction of this model is that the lensing galaxy, whose po¬ 
sition could be revealed by observing below the Lyman limit 
for this system ( ^ 4400A), hence removing the glare from 
the quasars, should be offset ~ 0.5arcsec from the quasar 
image, rather than lying behind the quasar images. 

Currently, our CO images of APM08279+5255 are of 
limited signal-to-noise. However, with further integration a 
detailed map of the CO image can be made. As this region 
will be free from the effects of microlensing, and as its ex¬ 
tended nature provides many more constraints (Kochanek, 
Keeton & McLoed 2001), such imaging has the poten¬ 
tial to provide a more accurate model of the lensing in 
APM08279+5255 than from the quasar images. 
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